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, Abstract. We present an ab initio calculation of the quark mass dependence of the ground state energies of 

Cs| ■ ■'He, *Be and ^^C, and of the energy of the Hoyle state in ^^C. These investigations are performed within 

the framework of lattice chiral Effective Field Theory. We address the sensitivity of the production rate of 
carbon and oxygen in red giant stars to the fundamental constants of nature by considering the impact of 
variations in the light quark masses and the electromagnetic fine-structure constant on the reaction rate 
of the triple-alpha process. As carbon and oxygen are essential to life as we know it, we also discuss the 
, implications of our findings for an anthropic view of the Universe. We find strong evidence that the physics 

^ ■ of the triple-alpha process is driven by alpha clustering, and that shifts in the fundamental parameters 

^ ' at the ~ 2 — 3% level are unlikely to be detrimental to the development of life. Tolerance against much 

larger changes cannot be ruled out at present, given the relatively limited knowledge of the quark mass 
dependence of the two-nucleon S'-wave scattering parameters. Lattice QCD is expected to provide refined 
^ estimates of the scattering parameters in the future. 

l/-^ ■ PACS. 21.10.Dr - 21.30.-x - 21.45.-v - 21.60.De - 26.20.Fj 
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' 1 Introduction while the total and radiative widths are r^^^. = 8.3(1.0) eV 

Q ,. and = 3.7(5) meV. 

[ While the Hoyle state dramatically increases the re- 
The "Hoyle state" is an excited state of ^^C with quan- action rate of the triple- alpha process, the resulting en- 
tum numbers = 0+ just above the ^Be-a threshold, hancement is also extremely sensitive to the exact value 
The existence of such a state was predicted by Hoyle in of e, which is therefore the principal control parameter 
order to explain the observed abundance of ^^C, which of this reaction. As the Hoyle state is crucial to the for- 
I is produced during helium burning in red giant stars via mation of elements essential to life as we know it, this 
■ ■ ■ the so-called triple-alpha process. In this two-step process, state has been nicknamed the "level of life" [1] Thus, the 
two ^He nuclei first combine to form the unstable, but rel- Hoyle state is often viewed as a prime manifestation of 
atively long-lived (under the conditions prevalent in the the anthropic principle, which states that the observable 
cores of red giant stars) ^Be nucleus. This resonance must values of the fundamental physical and cosmological pa- 
then combine with a third alpha particle in order to gen- rameters are restricted by the requirement that life can 
erate carbon. form to determine them, and that the Universe be old 
However, the fact that this process cannot by itself ex- enough for that to occur ^ . See, however, Ref. [7] for a 
plain the observed abundance of C in the Universe moti- thorough historical discussion of the Hoyle state m view of 
vated Hoyle in 1954 to propose the existence of an excited anthropic principle. We remark that in the context of 
0+ state of i^C, just above the ^Be-a threshold. Such cosmology and string theory, the anthropic principle and 
a resonant enhancement could then provide a sufficiently ^^s consequences have had a significant influence (see e.g. 
high rate of production of ^^C to account for the observed R-cfs. [HUH]). 

abundance [1]. Soon afterwards, the predicted state was The impact of changes in the energy of the Hoyle state 

detected at Caltech [HIS], and the modern value for its on the synthesis of ^^C and ^^O in red giant stars has 

energy is e = 379.47(18) keV relative to the 3a threshold, been investigated in several numerical studies that make 
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use of highly sophisticated stellar evolution models. Livio 
et al. |10| modified the value of e by hand and studied 
the triple-alpha process in the core and shell helium burn- 
ing up to the asymptotic giant branch stage in the stellar 
evolution. These calculations have been refined by Ober- 
hummer et al, who concluded that the production of ei- 
ther ^^C or ^^O becomes strongly suppressed for changes 
larger than 6{e) ~ ±100 keV in the position of the Hoyle 
state [TT] . In essence, if e is lowered too much, the triple- 
alpha process ignites at a significantly lower stellar core 
temperature, and hence not much energy is available for 
the process ^^C -I- ''He — >■ ^^O -I- 7. Conversely, if e is raised 
too much, the triple-alpha process ignites at a much higher 
core temperature, and hence most of the ^^C formed is im- 
mediately converted into ^^O and ^°Ne already before the 
conclusion of core He burning. However, since a ~ 100 keV 
change in e could still be tolerated, which is a 25% mod- 
ification, the degree of fine-tuning was revealed to not be 
as severe as was first believed [TJ ■ 

In addition to these ad hoc changes in e, a more mi- 
croscopic calculation was performed by Oberhummer et 
al. [T21IT3] in terms of a nuclear cluster model based on 
a simple nucleon-nucleon (NN) interaction with inclusion 
of electromagnetic (EM) effects. This NN interaction was 
formulated in terms of one strength parameter adjusted 
to give a fair description of a-a scattering and the spec- 
trum of ^^C. By modifying this coupling strength and the 
EM fine structure constant a^^^, the effect on the stellar 
production of ^^C and ^^O was analyzed. Within such a 
model, an adequate amount of ^^C and ^^O was produced 
within a relatively narrow window of ~ 0.5% around the 
observed strong force and of ~ 4% around the observed 
strength of the EM interaction. For larger changes, the 
stellar production of carbon and/or oxygen was found to 
be reduced by several orders of magnitude. 

However, the translation of the findings of Ref. [T5] 
into anthropic constraints on fundamental parameters re- 
mains problematic, as the employed model of the strong 
force is not readily connected to the fundamental the- 
ory of the strong interactions. Quantum Chromodynamics 
(QCD) and its fundamental parameters, the light quark 
masses. In this study, we address this question by means of 
an ab initio calculation of the sensitivity of e to changes in 
the light quark masses and the EM fine structure constant 
a^jjj. For this purpose, we carry out large-scale numerical 
lattice calculations for the energies and energy differences 
relevant to the triple-alpha process within the framework 
of chiral Effective Field Theory (EFT). The discretized 
(lattice) version of chiral EFT was formulated in Ref. [TS] 
(see Ref. [TB] for a recent review). We have successfully 
applied this novel approach to the spectra and properties 
of light nuclei [TTlfTSlfT^ I^ . to dilute neutron matter [H], 
and to the structure of the Hoyle state [22] . A brief sum- 
mary of the results reported here has appeared in Ref. j23j . 

This paper is organized as follows: In Section [21 we 
briefly describe the theoretical framework of the present 
calculation. The quark mass dependence of the nuclear 
force is described in Section [3] within the framework of 
chiral EFT, in combination with lattice QCD calculations. 



The analysis of the NN system is carried out in Section [JJ 
while Section [5| deals with the Auxiliary Field Quantum 
Monte Carlo (AFQMC) calculation of the energies of the 
^He, ®Be and ^^C ground states as well as the Hoyle state, 
including the energy differences relevant for the triple- 
alpha process. Section [S] provides an estimation of the 
neglected higher-order effects, while the observed corre- 
lations between the relevant energies and energy differ- 
ences are described in Section [7| Finally, in Section [5] the 
implications of our findings for the reaction rate of the 
triple-alpha process arc discussed, and Section [S] contains 
a concluding summary. 

2 Theoretical framework 

As discussed above, the triple-alpha reaction proceeds in 
two steps. The reaction rate of the first step "^He -I- ^He 
— )■ ®Be, where the unstable ®Be resonance is formed, is 
controlled by the energy difference AE^^, 

AE, ^E,- 2E„ (1) 

where we have introduced the notations E^ and E^ for the 
energies of the ground states of ^He and *Be, respectively. 
Also, we shall denote the ground state energy of ^^C by 
E12, and the energy of the Hoyle state by Ei2- The second 
step in the triple-alpha process, ^Be -I- ''He ^^C + 7, 
depends crucially on the energy difference AE/^, 

AE,, = EI2 -E,-E,, (2) 

such that the reaction rate for fusion of three a particles 
to '^^C via the ground state of ^Be and the Hoyle state is 
given by [13| 

where is the number density of a particles in a stellar 
plasma at temperature T. The energy difference e is given 

by 

e = AE, + AE,, = EI2 - SE,, (4) 

which clearly is the dominant control parameter of Eq. ([3|) 
in comparison with the linear dependence on the radiative 
width of the Hoyle state. The latter can therefore be 
neglected when the impact of small changes in e on r^^ 
are considered. 

The main question wc shall address here is the sensitiv- 
ity of e to variations in the light quark masses and the EM 
fine structure constant, with the objective of translating 
the bounds on e found in the stellar model calculations of 
Ref. [T^IM] to constraints on these fundamental parame- 
ters. For this purpose, we shall only consider the average 
light quark mass = [m^ + to^)/2, as the effects of 
strong isospin violation due to m„ 7^ are greatly sup- 
pressed for the processes of relevance to the present work. 
To this end, we shall calculate how E^, Eg^, E12 and E12 
depend on and aem, after which we may use this in- 
formation to determine the corresponding dependences of 
AEi^, AE,^ and e. 



Evgeny Epelbaum et al.: Dependence of the triple-alpha process on the fundamental constants of nature 



3 



Our analysis is carried out within the framework of chi- 
ral nuclear EFT, introduced by Weinberg [251 as a system- 
atic tool to explore the consequences of spontaneous and 
explicit chiral symmetry breaking of QCD in a rigorous 
manner. This approach relies on the most general effective 
Lagrangian for pions and nucleons constructed in harmony 
with the symmetries of QCD. The pions arc identified with 
the pseudo-Goldstone bosons of the spontaneously broken 
chiral symmetry of QCD, which strongly constrains their 
interactions. The small (but non-vanishing) pion mass is 
a result of the explicit breaking of chiral symmetry in 
QCD by the quark masses. In particular, one finds ^ 
{m^ + rrij), so that any dependence on the average light 
quark mass can be translated into a corresponding 
dependence on the pion mass M^. 

Chiral nuclear EFT is based on an order-by-order ex- 
pansion of the nuclear potential. In this scheme, two-, 
three- and four-nucleon forces arise naturally, and their 
observed hierarchy of importance is also explained. The 
nuclear forces have been worked out to high precision and 
applied successfully in few-nucleon systems to the bind- 
ing energies, structure, and reactions (see Refs. pSl[?7] for 
recent reviews). We have recently developed a discretized 
version of chiral EFT which allows one to compute the 
correlation function 



ZA{t) = {M^A<^M-tH)WA). 



(5) 



for A nucleons in Euclidean space-time using Monte Carlo 
sampling. Here, ^P'^ denotes the Slater determinant for A 
non-interacting nucleons, and H is the nuclear Hamilto- 
nian calculated in chiral EFT and expressed in terms of 
the lattice (discretized) variables. The correlation function 
Z^{t) can be efficiently calculated within the AFQMC 
framework, where terms in the lattice action quartic in 
the nucleon fields are re-exprcssed as interactions of a sin- 
gle nucleon with auxiliary fields by means of a Hubbard- 
Stratonovich transformation. Once Z j^(t) has been calcu- 
lated, the ground-state energy is obtained from the 
large-t limit of Zj^^ 



E, = -lu.'-^ 

^ t^oo dt 



(6) 



where cr.j and refer to the Pauli spin and isospin ma- 
trices of nucleon i, respectively. The strength g^j^ of the 
one-pion exchange potential is defined in terms of the nu- 
cleon axial- vector coupling and the pion decay constant 
F^as 

5 = ^ 

y-irN — 2p ' 



with t the Euclidean time. Additionally, we have devel- 
oped a multi-channel projection Monte Carlo method, which 
allows us to study excited states by computing the corre- 
lation matrix for a set of A-nucleon states with appro- 
priately chosen quantum numbers. 



z^i{t) = {n\eM-tH)\n)- 



(7) 



The AFQMC results reported here correspond to an 
improved leading-order (LO) action, based on the NN am- 
plitude 



-^LO = C'o /(g) \ 



+ C,/(g)g + i<T,.^,. 



3 1 
1 1 



T • T 

4 4' 



— — — X • T 

4 4 » J 



(8) 



(9) 



while Cq and Cj are low-energy constants (LECs), ad- 
justed to reproduce the NN phase shifts in the ^S^ and 
^S^ partial waves. The smearing function /(q) is chosen 
to give the (approximately) correct effective ranges for 
the two S-wave NN channels (see Ref. [TH] and references 
therein for more details, along with a description of the 
discretized form of the improved LO action) . Higher-order 
corrections are taken into account in perturbation theory. 
In particular, expectation values of a given operator O are 
obtained as 

Zl{t) = (S'aI cxp(-tF/2) O cxp(-<F/2)|if-^), (10) 

which accounts for all contributions to the nuclear Hamil- 
tonian up to next-to-next-to-leading order (N^LO) in the 
chiral expansion, including the Coulomb interaction and 
the threc-nucleon forces. 

Our recent AFQMC calculations using this framework 
are reported in Refs. pn[r8..19..20.,22] . In particular, re- 
sults at N^LO for nuclei with A = 3,4, 6 and 12 can be 
found in Ref. [TH]. These calculations employ a periodic 
cubic lattice with a lattice spacing of a = 1.97 fm and 
a length of L = 11.82 fm. In the (discretized) Euclidean 
time direction, we use a step size of a^ = 1.32 fm, and 
perform calculations for propagation times, i.e. the ex- 
tent of the time direction, = N^a^ = 4 ... 20 fm, such 
that the limit — >■ oo is taken by extrapolation. Given 
the relatively coarse lattice spacing employed in our cal- 
culations, the two-pion exchange NN potential that starts 
contributing at next-to-leading order (NLO) can be well 
represented by contact interactions |19j . 

We now turn to the m^-dependence or, equivalently, 
the M^-dependence of the energies E^ . We use the nota- 
tion E^ when referring to either the energies E^^ E^, E12 
and E^2: or to the energy differences AEf^, "^^^h ^^'^ ^■ 
In this work, we shall restrict the values of to the 
vicinity of the physical pion mass, roughly speaking to 
\SM^/M^\ < 10%. It is then sufficient to consider the lin- 
ear variation of E^, giving 



SE, 



dE, 



dE, 



MS 



(11) 



where we have allowed for independent variations of 
and aom- Our objective is then to compute the partial 
derivatives in Eq. ([TT|) using AFQMC. Clearly, such a cal- 
culation relies on knowledge of the M^-dependence of the 
nuclear Hamiltonian, which is discussed in Section [3l We 
also note that a useful way to express the sensitivity of 
a given observable X to a, parameter y is given by the 
dimensionless "if-factors" 



y_ dX_ 
X dy 



(12) 
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where we use the superscript i = {q, tt, a} for the set of 
observables y = {m^, M^, ctcm}- example, we can 

obtain Kj^ (i.e. the sensitivity of X to changes in m^) in 
terms of by means of the relation 



0.494 



+0.009 
-0.013 



(13) 



from 



where we shall adopt the value K\i 

Ref. ng. 

In addition to shifts in m^, we shall also consider the 
effects of shifts in a^^. The treatment of the Coulomb 
interaction in our AFQMC framework is described in de- 
tail in Ref. [T2]. The main difference between the con- 
tinuum and latticc formulations is that the discretizcd 
form of the long-range Coulomb force between two pro- 
tons becomes singular if the protons occupy the same lat- 
tice site. We therefore employ a regularized version of the 
discretized Coulomb interaction, where the potential en- 
ergy of two protons on the same lattice site is set to the 
continuum value corresponding to a separation of half a 
lattice spacing. The effects of this regularization are com- 
pensated for by a derivative-less proton-proton contact op- 
erator, which also receives contributions from the strong 
and short-range EM isospin-breaking effects. The associ- 
ated coefficient c^^ is determined from the proton-proton 
phase shifts. 

The sensitivity of the energies to variations in 
can be obtained by computing the shifts AE^{a^^ and 
AE^{Cpp). The former is due to the long-range Coulomb 
interaction on the latticc, and the latter arises from the 
part of the proton-proton contact operator oc Cpp. Specif- 
ically, we define 



3 Pion mass dependence of the nuclear 
Hamiltonian 

The Af^-depcndence of the energies E^ is generated by the 
M^-dcpendence of the nucleon mass in the kinetic en- 
ergy term in the nuclear Hamiltonian, as well as by the 
M^-dependence of the nuclear potentials. In the present 
analysis, we will not take into account the sources of M^- 
dependence generated by contributions beyond Aj^q in 
Eq. dS]). Instead, we shall estimate the neglected higher- 
order contributions to the energy shifts in Section [S) On 
the other hand, in order to obtain the most accurate de- 
scription possible of the M^-dependence of .4^0 ^^'^ ^'^ 
ensure model independence, we shall go beyond the strict 
chiral expansion of the terms entering Eq. ([8]), and make 
use of the available lattice QCD data whenever possible. 

For the M^-depcndence of the nucleon mass, we define 
the quantity 



9A/ 



(17) 



which has been analyzed extensively in the literature by 
combining Chiral Perturbation Theory (ChPT) with lat- 
tice QCD data, see e.g. Ref. ^^M- At 0{p^) in ChPT, 
"%(^'^7r) is given by 



(18) 



Q,^{E^) = AE,{r 



(14) 



where ttiq denotes the value of in the chiral limit. The 
value ~ 0.57 was obtained in Ref. [25]: which corre- 
sponds to the 0{p'^) heavy-baryon (HB) ChPT result with 
Ci = -0.81 GeV^^ Alternatively, one may determine x^ 
from the pion- nucleon sigma term cr^^ = M^dniN /dM^ 
by means of the Feynman-Hellmann theorem. From the 
results of Ref. [33] , one finds 



where the coefficient 



denotes the relative strength 



of the pp contact interaction caused by the regulariza- 
tion of the Coulomb force. The determination of Xpp will 
be described in Section [5] Given that the energy shifts 
AE^{a^^) and AE.^{Cpp) are relatively small, we may ap- 
proximate 



ph 



(15) 



with ~ 1/137. While we shall mainly study the ex- 
plicit dependence on induced by the Coulomb inter- 
action, it is worth noting that K'^ may receive additional 
contributions from the corresponding shifts in the effec- 
tive hadronic Lagrangian, for example from the EM shift 
of M^. Schematically, this may be expressed as 



i\x — 



Q,^{X) , AMr dX 



X 



X aM^ 



(16) 



where the sizes of such additional contributions to the EM 
il'-factors are Z\M™/M^ ~ 4% compared to the domi- 
nant strong contributions, and are thus not considered in 
the present analysis. 



a^^ = U.Qtli MeV 



while those of Ref. [SHj yield 



CT^^ = 59 ± 7 McV 



•^1 = O-^^^^-o^ol (19) 



= 0.87 ±0.10, (20) 



and we also note that Xi ~ 0.73 has been obtained by 
Procura et al. from fits of a modified 0{p'^) ChPT formula 
to lattice QCD data [52]. In our analysis, we adopt the 
conservative estimate 



0.57... 0.97, 



(21) 



based on the variation in the results quoted above. 

We now turn to the M^-dependence of the nuclear 
force. The most obvious source of M^-dependence is the 
static pion propagator in Eq. ([5]). In addition to this ex- 
plicit dependence on M^, we also take into account the 
implicit A/^-dependence of the coupling constant g^^ of 
the one-pion exchange (OPE) potential by defining 



N 



ME 



1 

2f: 



dg^ 



MS 



_9a_ dK 
2Fi dAE 



MS 



(22) 
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where both contributions to X2 have been studied exten- 
sively by means of ChPT and lattice QCD (see Refs. [28l 
[5D1I5T] and references therein). As for in Eq. p^ . we 
obtain 



Ml 



Wtt^F'^ 



(23) 



where we use ~ 4.3 from Ref. 



which is consis- 
see e.g. Ref. [55] . 



tent with modern lattice determinations, 
Ti = 4.67(3)(10), and F ~ 86.2 MeV is the value of F^ in 
the chiral limit. As an alternative to the sub- leading or- 
der ChPT result, one may use the determination of K'^ 

reported in Ref. ^ 0.048 ± 0.012, which is based 

on a combined analysis in ChPT and lattice QCD. This 
gives 



dF^ 



MS 



0.066, 



(24) 



using the central value K\j = 0.494, as discussed in Sec- 
tion n 

While the chiral expansion of F^ shows good conver- 
gence, the equivalent expression for is known to con- 
verge much slower. In particular, the 0{p^) HB ChPT re- 
sult for shows a very strong M^-dependence near the 
physical point (sec Refs. p8l[34] and references therein). 
On the other hand, lattice QCD calculations indicate that 

as a function of is remarkably flat. In principle, 
such flat behavior can be accommodated at the two-loop 
level in ChPT [55]. By assuming that 5^ is constant in 
the vicinity of the physical point, the value X2 — —0.049 
in lattice units (l.u.) is obtained from Eq. ([5^ . This should 
be compared with the central value reported in Ref. ^5] . 
X2 = —0.024 l.u., which is obtained from the incomplete 
0{p^) HB result for g^, constrained by the available lat- 
tice QCD result at ~ 350 MeV. The relatively large 
uncertainty in the chiral extrapolation of has a signifi- 
cant impact on the allowed values of X2 ■ The largest source 
of uncertainty in Ref. |28| is the poorly known low-energy 
constant d^Q. The present empirical constraints from the 
reaction ttA^ — !■ ttttN yield a relatively large range of 



^16 



-0.91. 



-2.61 GeV , which in turn gives 



X2 = -0.056. 



.0.008 l.u. 



(25) 



for the range of uncertainty in X2, which we shall adopt 
in the present analysis. 

The short-range part of the nuclear force also depends 
on . This dependence is more difficult to control within 
chiral EFT, see Ref. [55] for an extended discussion. Since 
we aim at a model-independent determination of the M^- 
dependence of the nuclear energies E^, we refrain from 
a chiral expansion of the short-range part of the nuclear 
force. Rather, we shall parametrize the M^-dcpendcncc of 
the LO contact interactions, i.e. of the coefficients Cg and 
Ci in Eq. ([5]). This can be performed in terms of the slope 
of the inverse NN S'~-wave scattering lengths aj^ and a^^. 



where we have introduced the subscripts s and t for the 
spin-0 Q'Sq) and spin-1 QSi) NN partial waves. For the 
purpose of our analysis, and are regarded as input 
parameters, and we shall express all our results in terms 
of these. We shall return to the determination of A^ and 
A^ in Section [5] 

Given the sources of M^-depcndence discussed so far, 
we may express the dependence of the energies E^ on 
as 

E, = ii;,(M,,m^(Mj,g,^(MJ,Co(A4),C,(Mj), 

(27) 

where refers to the explicit M^-dependence from the 
pion propagator in the OPE contribution. In order to as- 
sess the sensitivity of the triple-alpha process (and of the 
various energy levels involved in that process) to shifts in 
M^, we will compute quantities of the form dE^/dM^ at 
the physical point. Given Eq. ([57)1 . we find 



dE, 



OAF 



dE, 



dE, 



dM. 
dE, 



Ml 



^ dm 



N 

dE., 



dE. 



N 



dCj 



where we have introduced the definitions 



9Cn 



OAF 



dCj 



MS 



dhF 



(28) 



(29) 



MS 



for the short-range components of the LO amplitude. Our 
method for re-expressing the scheme-dependent parame- 
ters and X4 in terms of A^ and A^ is explained in Sec- 
tion H] The AFQMC calculation of the partial derivatives 
in Eq. (j28[) is detailed in Section [5] 

4 Short-range contributions to the nuclear 
force 

With the exception of the OPE contribution, much of the 
Af^- dependence of the nuclear Hamiltonian is implicit and 
thus controlled by the coefficients Xj, which describe how 
mjY,5^jY,CQ and Cj depend on M^. As discussed in Sec- 
tion [3] Xi and X2 are fairly well constrained by ChPT cal- 
culations in combination with lattice QCD data, and pro- 
vide external input to our analysis according to Eqs. (|2ip 
and (|25p . Contrary to x^ and X2, the coefficients X3 and 
x^ are strongly scheme-dependent. We shall therefore ex- 
press these in terms of the scheme-independent quantities 
A^ and A^ . This substitution is achieved by an analysis of 
the two-nucleon scattering problem on a periodic lattice 
of length L. Once the dependence of the S'-wave scatter- 
ing lengths in the spin-singlet and triplet channels on 
is known, x^ and x^ can be straightforwardly obtained. 
Specifically, we make use of the finite volume formula due 
to Liischer 



A^ = 



da. 



dAF 



M^ 



MS 



(26) pcotS^^Siv)^--, 
ttL a 



(30) 
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which we can differentiate with respect to M^. This yields 



(31) 



where the function S{'q) and its derivative may be com- 
puted numerically using the expansions 



(32) 



with the coefficients S'„ given by 

5o = -8.913631, 
^2 = 8.401924, 
^4 = 6.426119, 
^6 = 6.098184, 



S*! = 16.532288, 



6.945808, 
6.202149, 



(33) 



= 6.048263. 



If we denote the two-nucleon energy in the singlet channel 
by Eg and that in the triplet channel by E^, we find 



2n r-*ail/. 



(34) 



and by introducing the notation of Eq. (|26p , we obtain the 
relations 




+ (Xg - 3X4) 



where we have defined 



- 4^3 ^^«,t''' 



dE^ 



<ls,t 



5Cn 



(37) 



Table 1. LO energies (in MeV) of the spin-singlet [Eg) and 
spin-triplet (iSj) two-nucleon states used in the Liischer analy- 
sis, along with the required partial derivatives. The results were 
obtained by numerical solution of the Schrodinger equation in 
a cubic box of size N = 2A (second column) and A'^ = 32 (third 
column). denotes the energy of the deuteron. Quantities 
labeled "[l.u.]" are given in units of the inverse (spatial) lattice 
spacing. All derivatives are evaluated at the physical point. 





L = 47.36 fm 


L = 63.14 fm 


E,{hO) [MeV] 


-0.0440956 


-0.0218593 


£;t(LO) [MeV] 


0.0369820 


0.0142463 


E^{hO) [MeV] 


-2.2244401 


-2.2443719 


BEJdM^ 


-2.0261727x10"'' 


-1.1927443x10"'' 


dEJdM^ 


-5.5518808x10"^ 


-1.9921726x10"^ 


dEg/drriff 


-2.0081283 xlO""* 


-1.2224536x10"* 


dEf/drrif^ 


-1.6013626x10""* 


-5.8480992x10"^ 


dEJdg^N[\.u.\ 


9.9946261x10"* 


5.8797858x10"* 


dEtldg^N[^.\i.\ 


3.5402207x10"'' 


1.2691192x10"* 


qa [1-u.] 


0.00379650 


0.00223108 


qt [l.u.] 


0.00165886 


5.9467427x10"* 



parameters contains sizable uncertainties. Elimination of 
and in favor of and Af gives 



dE, 



Q 



MC 



MC 



Qt ^2) + ^Mc(-^l' ^2)' (38) 



which is equivalent to Eq. ([28]) . Here, the notation "MC" 
indicates which quantities incorporate information from 
the AFQMC framework. In this parametrization, the in- 
dividual terms may be obtained using the numbers in Ta- 
ble m together with the AFQMC results of Section [5l In 
Eq. ([55]) . we have introduced the quantities 



qMC ^ _^ 

J_ dE^ 



1 <9K 



pph 4q, dCj 

_±dE^ 

(^ph Aq^ dCj 



(39) 
(40) 



Our results for the energies i?^ ^ and the corresponding 
partial derivatives, including the factors ^, are summa- 
rized in Table [TJ These are computed by exact numerical 
solution of the two-nucleon problem on a spatial lattice. 
As the objective is to take the box size N (with L = Na) 
large enough to make finite volume effects negligible, two 
different box sizes {N = 24 and N = 32) have been consid- 
ered in order to determine the magnitude of residual finite 
volume effects on the Liischer analysis. Given Eqs. (|35p 
and and the results in Table [U we are in the position 
to compute x^ and x^ for use in Eq. 



However, it is 

also instructive to eliminate x^ and x^ analytically, which 
provides an alternative to the parametrization of Eq. psp . 
This is particularly useful when we express our final results 
in terms of A^ and A^, as the current knowledge of these 



for which the error is given entirely by the statistical un- 
certainty of the AFQMC calculation. We also have 



Rf{xi, X2) 




(41) 



(42) 
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where, in contrast, the dominant sources of uncertainty 
come from Xi and X2- Finally, we have 



dm 



N 



dE^ 



N 



(43) 

which combines the AFQMC results for the OPE and ki- 
netic energy contributions. The error of Ry^Q receives con- 
tributions from the statistical AFQMC error as well as 
from Xi and X2- 

While the most convenient way to obtain our final re- 
sults is by means of Eq. ([55]) , we may also use the values 
given in Table [T] to solve Eqs. ([55]) and ([M]) for Xg and x^, 
and to express these as functions of ^ and Xi2- This 
also allows us to illustrate the size of the finite volume 
effects in the Liischer analysis. We obtain the relations 



4.8394 X 10^2 + 6.7146 x IQ'^ x^ 



(44) 



- 0.25080 - 0.37540^, - 0.20377 A^, 
x^ = 4.9754 X IQ-^ - 1.8813 x 10"^ x^ (45) 
1.2462 X 10-^X2 -0.12513^, + 0.20377 Aj, 



for the smaller = 24 lattice {L = 47.36 fm), and 

= 4.8470 X 10-2 + 6.7127 x 10-^ x^ (46) 

- 0.25101 X2 - 0.37652 - 0.20467 A^, 

Xi = 4.9901 X 10-^ - 1.8998 x 10"^ x^ (47) 

- 1.2532 X 10-2x2 - 0.12551 + 0.20467 At, 

for the larger A^ = 32 lattice [L = 63.14 fm). In the 
above equations, the dimensionful quantities X2, x^ and 
should be taken in units of the corresponding powers 
of the inverse lattice spacing. We note that the results for 
A' = 24 and A = 32 are practically indistinguishable. In 
Section [5l we shall make use of the A^ = 32 lattice when 
presenting our final results. 



5 Auxiliary Field Quantum Monte Carlo 
results 

We now turn to the AFQMC calculation of the shifts in 
the E^. Our Monte Carlo simulations arc performed for 
a single value of M^, equal to the neutral pion mass, 
with isospin symmetry breaking treated as a perturba- 
tion. The partial derivatives dE^/dC^ and dEjdCj are 
obtained by computing the matrix elements of the as- 
sociated operators according to Eq. (flOl) . On the other 
hand, the partial derivatives dE^/dM^ with respect to 
the pion mass in the OPE term are computed by evaluat- 
ing in perturbation theory the energy shift AE^ induced 
by the substitution H{M^) H{M^ + AM^) in the nu- 
clear Hamiltonian. Here, the masses of both the neutral 
and charged pions in the OPE term have been shifted by 
AM^ = 4.59 MeV, which equals the empirical mass dif- 
ference between the neutral and charged pions. The corre- 
sponding partial derivatives dE^/dM^ that enter Eqs. (j28|l 



and (|43p are then given by 



dE^ 



MS 



^E,{AM^) 
AAf 



(48) 



which we find accurate to within the statistical error of 
the AFQMC calculation. We shall also briefly consider 
the closely related energy shift AE^{AAd^), given by the 
substitution H{M^) H{M^ + AM^) in the nuclear 
Hamiltonian. In this case, only the masses of the charged 
pions are shifted to their physical value. Finally, the par- 
tial derivatives dE^/dg^j^ and dE^/drhj^ are obtained as 
a finite difference, by defining the quantities g^jy ± Ag^j^ 
and mj^dzAnij^, followed by computation of the resulting 
shifts of the E^ in perturbation theory. 

All AFQMC results presented here have been extrapo- 
lated to infinite Euclidean time {N^ = 00). Such an extrap- 
olation is necessary, as increasing the number of Euclidean 
time steps beyond = 14 for '^Hc and Nf. = 12 for the 
heavier nuclei becomes impractical due to the worsening 
sign problem. An accurate extrapolation is necessary for 
reliable conclusions. As an example, on an A^ = 6 lattice 
the calculated *He binding energy at = 14 still devi- 
ates from the extrapolated value at the ~ 10% level. We 
observe a significantly faster rate of convergence with 
for the heavier nuclei, however, we also note that this de- 
pends strongly on the choice of trial wave function. The 
extrapolation of the LO energy E^ is performed using 



E,{N,)^E^ioo)+CEexp{- 



(49) 



and all matrix elements computed in perturbation theory 
(energy shifts and partial derivatives, collectively labeled 
Aj) are extrapolated with 



^.(^J = ^.M+c,exp 



2r 



(50) 



using a common exponent r. The extrapolation is per- 
formed by means of a simultaneous chi-square minimiza- 
tion for all quantities, such that the fitted parameters arc 
i?^(oo), Aj(oo), c^, Cj and r. We find that the stability of 
such an extrapolation requires that the data for the matrix 
elements be excluded for < 6, except for ^^C. 

As the extrapolation to A'j = 00 is crucial for the ac- 
curacy of our conclusions, it is highly desirable to validate 
our method against known results. Such a comparison is 
given for the deuteron in Table where the extrapo- 
lated AFQMC results (see Fig. [IJ are compared with a 
direct numerical solution of the Schrodinger equation for 
an identical Hamiltonian. The AFQMC test results for 
the deuteron were obtained in a relatively small box size 
of i = 5.92 fm (A^ = 3). The comparison in Tabic [5] shows 
excellent agreement between the extrapolated and exact 
results, which gives confidence that the extrapolation pro- 
cedure employed in our analysis is indeed reliable. We may 
then proceed with the AFQMC calculation for the nuclei 
relevant to the triple-alpha process. The results for ^He 
and ®Be are summarized in Table [3J and the individual 
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Table 2. Validation of the extrapolation A'^t — i- cxj to infinite 
Euclidean time. The LO deuteron energy and the corre- 
sponding energy shifts and derivatives (at the physical point) 
are computed using AFQMC and extrapolated ^ oo (sec- 
ond column) and compared with the values obtained from the 
numerical solution of the Schrodinger equation (third column) . 
The appropriate units are given for each quantity, with "[l.u.]" 
indicating units of the inverse (spatial) lattice spacing. Paren- 
theses indicate one-standard-deviation errors. 





(MC-fex) 


"^H (exact) 


Ea{LO) [MeV] 
AE^iAM^) [MeV] 
AE^{AM'^) [MeV] 
dEJdmj^ 
dEa/dg^M [l.u.] 
dEJdCo [l.u.] 
dEJdCr [l.u.] 


-9.070(12) 
-0.003548(12) 
-0.002372(8) 
-0.00382(2) 
0.01024(11) 
0.13897(15) 
-0.4171(4) 


-9.078 
-0.003569 
-0.002379 
-0.003809 
0.01017 
0.138867 
-0.41660 



Table 4. AFQMC results for the Hoyle state (second column) 
and the ground state (third column) of ^^C. The LO ener- 
gies Ei and the corresponding energy shifts (including the EM 
shifts) and derivatives (at the physical point) have been extrap- 
olated Nt ^ oo. The appropriate units are given for each quan- 
tity, with "[l.u.]" indicating units of the inverse (spatial) lat- 
tice spacing. Parentheses indicate one-standard-deviation er- 
rors. All derivatives are computed at the physical point. 





^^C* (MC-fex) 


'^C (MC+ex) 


E,{LO) [MeV] 


-89.8(13) 


-95.6(6) 


AE,{AM^) [MeV] 


-0.802(2) 


-0.778(4) 


AE.ic^^) [MeV] 


2.032(10) 


1.95(2) 


AE,{a,^) [MeV] 


5.54(2) 


5.67(2) 


dEi/dniff 


-0.403(5) 


-0.395(5) 


dEJdg^N [1-U-] 


1.343(13) 


1.285(16) 


dEJdCo [l.u.] 


6.86(3) 


6.55(7) 


dEJdCj [l.u.] 


-7.92(3) 


-7.54(6) 



Table 3. AFQMC results for ''He and **Be. The LO ener- 
gies Ei and the corresponding energy shifts (including the EM 
shifts) and derivatives (at the physical point) have been extrap- 
olated ^ oo. The appropriate units are given for each quan- 
tity, with "[l.u.]" indicating units of the inverse (spatial) lat- 
tice spacing. Parentheses indicate one-standard-deviation er- 
rors. All derivatives are computed at the physical point. 





*He (MC-l-ex) 


''Be (MC4-ex) 


Ei(LO) [MeV] 


-28.89(11) 


-57.2(5) 


AEi{AM^) [MeV] 


-0.2290(17) 


-0.477(5) 


AEiiCpp) [MeV] 


0.433(3) 


1.02(3) 


AEiia,^) [MeV] 


0.613(2) 


2.35(2) 


dEi/drrif^ 


-0.0750(7) 


-0.187(6) 


dE,/dg^N [l.u.] 


0.337(3) 


0.746(12) 


dEJdCo [l.u.] 


1.527(12) 


3.52(8) 


dEJdCj [l.u.] 


-1.881(17) 


-4.22(7) 



MC data points along with the extrapolation are shown in 
Figs. [5] and O The results for the ground and Hoyle states 
of '^C are similarly given in Table 2] and Fig. These 
AFQMC calculations were performed on a L = 11.84 fm 
{N = 6) lattice, which is large enough to render residual 
finite volume effects smaller than the expected error from 
truncation of the chiral expansion of the NN interaction 
at N^LO. 

We may now combine our AFQMC results with the 
two-nucleon scattering analysis, in order to obtain pre- 
dictions for the M^-dependence of the various states fea- 
turing in the triple-alpha process. This can be performed 
straightforwardly by substituting the AFQMC results and 
the numbers from Table [1] into Eq. ((38)) . and by propagat- 
ing the various sources of error. As described in Section|31 
we adopt 



73+0-24 



— U.UZ4_Qg34 



Lu., 



(51) 



for the central values and uncertainties in Xi and X2, and 
express the results as a function of and A^. In this way, 
we obtain the following results for the M^-dependencc of 



the energy levels involved in the triple-alpha process, 



dE,, 



MS'' 



ME 



dEt, 



MS 



- 0.339(5) A, - 0.697(4) A^ 
4-0.0380(14)t°:°°8, 

- 0.794(32) A^ - 1.584(23) A^ 
+ 0.089(9)±°:^il, 

- 1.52(3) ^ - 2.88(2) At 
+ 0.159(7)t°;°fs, 

- 1.588(11) - 3.025(8) At 
+ 0.178(4)t°;^^6^ 



(52) 



(53) 



(54) 



(55) 



where the error receives contributions both from the sta- 
tistical error of the AFQMC calculation (given in paren- 
theses) as well as from the uncertainties in x-^ and (ex- 
plicit positive and negative bounds given). It is notewor- 
thy that Xi and X2 only affect the constant terms in the 
above results, and therefore the (sizable) uncertainties in 
these coefficients have a relatively minor impact. We can 
also assess the sensitivity to small shifts in by comput- 
ing the "iC-factors" as defined in Eq. (fT^ . For this pur- 
pose, we take = 138.0 MeV as the isospin-averaged 
pion mass, and the empirical values E^^^ ~ —28.30 MeV, 
£:^'^P = -56.50 MeV, E^^^ = -92.16 MeV, and E*^''^ = 
-84.51 MeV for the E^. This yields 



1.652(25) A, 



1.94(8) A, 



- 3.401(21) At - 0.185(7)1° °^^ 



0.039' 

(56) 

-1-0.027 



K 
K 



3.87(6) At- 0.217(21)1° (57) 
2.27(4) A, + 4.32(3) A^ - 0.239(11)1^;^3^, (58) 



2.593(19) A, + 4.940(13) A^ - 0.291(7)l^;^i^3, 



(59) 
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where the same conventions for the errors have been ap- 
phed. Having calculated the shifts of the individual energy 
levels involved in the triple-alpha process, we may com- 
bine these and obtain similar predictions for the energy 
differences AEf^, AEf^ and in particular for e, which is the 
critical control parameter for the triple-alpha reaction rate 



'3a 



in Eq. ©. We find 



dAE, 



dAE, 



de 



MP 



MP 



-0.117(34) A^- 0.189(24) 

+ 0.013(9)™, 

- 0.455(35) - 0.744(24) A^ 
+ 0.051(10)™, 

- 0.572(19) A^- 0.933(15) 

+ 0.064(6)™, 



and 

K 



--175(51)^- 284(36) A, 
= -217(16)^- 355(12) A, 
A7 ^ -208(7) A^ - 339(5) A^ + 23(2) 



19(13)t^:g, 



25(5)^1°, 



+3.7 
-3.4' 



(60) 
(61) 

(62) 

(63) 
(64) 
(65) 



where we have used the empirical values AE'^^^ = 92 keV, 
289 keV, and e = 380 kcV. 
We now turn our attention to the AFQMC results for 
the EM shifts. As explained in Section [51 our objective is 
to compute the energy shifts Qcmi^i) defined in Eq. (fH]) . 
This involves the unknown parameter Xpp, which deter- 
mines the relative strength of the proton-proton contact 
interaction that emerges from the lattice regularization of 
the long-range Coulomb force. We may fix by means of 
the known contribution of the Coulomb force to the bind- 
ing energy of ^Hc, which is Q^^iE^) = 0.78(3) MeV [S^. 
The quoted error refiects the model-dependence and is de- 
termined from the range of values corresponding to differ- 
ent phenomenological two- and three-nucleon potentials 
as well as nuclear forces derived in chiral EFT. Using the 
AFQMC results for AE^{a^„^) and AE^{Cpp) from Tabled 
we find 

Q^^{E^) = 0.433(3) MeV x x^^ + 0.613(2) MeV 

= 0.78(3) MeV ^ x^^- 0.39(5), (66) 

for ^Hc, which enables us to predict 

Q^^^iE^) = 1.02(3) MeV x Xpp + 2.35(2) MeV 

= 2.75(8) MeV, (67) 

for *Be. Further, using the value Xpp = 0.39(5) and the 
AFQMC results given in Table HI we predict 

<9cm(^^i*2) = 2.032(10) MeV x x^p + 5.54(2) MeV 

^ 6.33(6) MeV, (68) 



i-n included 
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— fit V*"^ 
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Fig. 1. AFQMC calculation of the deuteron, as a function of 
Euclidean time A^^. Results after the extrapolation A^f — >■ oo are 
given in Table [J] The results for Ej{LO) have been omitted, 
as they show no appreciable dependence on N^, and can thus 
be trivially extrapolated. 



for the Hoyle state, and 

(9cm(^i2) = 1-95(2) MeV x Xpp + 5.67(2) MeV 



6.43(6) MeV, 



(69) 



for the ground state of ^^C. We are now in the position to 
predict the EM shifts of the energy differences AEf^, AEf^ 
and e. This gives 



Q,^{AE^) - 2.80(10) MeV, 
Q,^{AE,) = 1.19(8) MeV, 
Q„Je) = 3.99(9) MeV. 



(70) 



We have also studied the dependence of the excitation 
energy AE^ of the Hoyle state on and Og^. Although 
AE^ is not needed for the calculation of the triple-alpha 
reaction rate, it provides an instructive reference point 
when discussing the sensitivity of various energy differ- 
ences to small changes in the fundamental constants. We 
compute AE^ from 



AE, EE Et, - E, 



(71) 
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Fig. 2. AFQMC calculation of ^He, as a function of Euclidean 
time steps Nt- Results after extrapolation N-t ^ oo are given 
in Table [3 The results of Ref. [TS] for E^^iLO) are included to 
highlight the improved statistics, and as a consistency check. 



for which wc find 
dAE, 



0.07(3)^,-0.14(2)^4 



MS 



0.019(9)Io 



0.004 



This yields 



Kae = -1-3(5) 4 - 2.6(4) + 0.34(15)i" "^ 



0.05 ' 



(72) 



(73) 



for the sensitivity to small changes in M^. The above 
result corresponds to the empirical value of AE°^p = 
7.65 MeV. Finally, we find 



Q,^{AE,) = 0.10(7) MeV, 



(74) 



Fig. 3. AFQMC calculation of *Be, as a function of Euclidean 
time steps Nt- Results after extrapolation N-t ^ oo are given 
in Table [3] The results of Ref. jS] for -E'g(LO) are included to 
highlight the improved statistics, and as a consistency check. 

for the EM shift in the excitation energy of the Hoyle 
state. 



6 Theoretical uncertainties and higher-order 
corrections 

In order to estimate the theoretical uncertainty of our re- 
sults, we shall first consider the effects of neglected higher- 
order terms in the chiral EFT expansion. To this end, we 
compute the Af^-dependence of the ^He binding energy 
induced by the explicitly A/j^-dependent part of the three- 
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nuclcon force (3NF), 



/iSNF 



9\ 



• X, (-4ciA/2 



8F4 (q2 + M2)(q2 + M2) 
9aD Qi -r, ■ 



permutations, 



Ml 



E 

— T • T ■ 

2 * ^ 



(75) 



which contributes at N'^LO. For the LECs Cj, we adopt the 
values determined from low-energy pion-nucleon scatter- 
ing data, namely = -0.81 GeV"\ Cg = -4.7 GeV"\ 
and C4 = 3.4 GeV"! [31131]. The LECs D and are fixed 
by means of the triton binding energy and the weak ax- 
ial vector current. For more details on the treatment of 
the 3NF in the nuclear lattice simulations, see Ref. [TJ- 
As a representative example, we calculate the sensitivity 
of to changes in entering the 3NF by performing 
AFQMC calculations using a slightly shifted pion mass 
in Eq. ([75]), namely = M^^ ± 20 MeV. The resulting 
shifts in E^ induced by the 3NF are 



4 IaA =A/P"+20 McV 



^^f ^Im =Mr'-2o McV = -3.204(16) McV, (76) 



= -3.034(12) MeV, 



which gives us the rough estimate 



0.004 



(77) 



This is an order of magnitude smaller than e.g. the LO 
contribution dE^/dM^\ j^,jph ~ —0.05. This observation 
suggests that the 3NF effects are indeed very much sup- 
pressed, as expected for a N^LO contribution. 

It is also instructive to compare our results with the 
ones of Ref. [51j . which were obtained within the frame- 
work of planless EFT. Not only does this provide a useful 
consistency check for our calculations, but it also allows 
us to estimate our theoretical uncertainty in a comple- 
mentary way. In particular, we can compare our result for 
dE^/dM^ in Eq. with the result for the ''He binding 
energy B4 {B^ = \B^\ = ~E^) given in Eqs. (1.5)-(1.7) of 
Ref. [S] 



0.037-^^^"^ 



0.74^^^^' 



(78) 



where B^ is the deuteron binding energy. Noting that 
d/dm^ oc d/dM^ and making use of the relation 
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Fig. 4. AFQMC calculation of ^^C, as a function of Euclidean 
time steps Nt. Results after extrapolation N-t ^ 00 are given 
in Table H The results of Ref. [H] for £T2(LO) are included to 
highlight the improved statistics, and as a consistency check. 

where it is still necessary to convert the dependence on 
B^ into a corresponding dependence on a^. To this end, 
we use the effective range approximation 



pcotd^ ~ + -p r^, 



(81) 



da. 



da. 



we can bring Eq. (|78p into the form 
dB^ „ da7^ 



-0.037 5,0 — 

or 



0.74 



B^dB^ 



(79) 



(80) 



to obtain 
dB, 



da;' 



dM^ [a, + ^a^{a^-2r,) - 2r,] dM^ 

daT^ 



: 0.164 



dAL 



(82) 
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where we have inserted the empirical values ~ 5.42 fm 



and 



1.75 fm. Finally, by eombining this expression 



with Eq. dini) and noting that = the pionless EFT 
result of Ref. [ST] is brought into the form 



0.126 A, -0.739 A 



E^m^af 



■A 



(83) 



which is in a reasonable agreement with Eq. (|52|) . By 
taking the central values of the coefficients Xi and X2 in 
Eq. (IS21), as well as = 0.29 and = -0.18 as found 
in the most recent chiral EFT calculations (see Section |S] 
for details), we obtain dE^/dM^ ~ 0.065. This should be 
compared with the pionless EFT result dEjdM^ ~ 0.096 
based on Eq. ([55)1 . We expect that the theoretical uncer- 
tainty of our calculation is much smaller than the differ- 
ence between these two numbers or. more generally, than 
the difference between Eqs. (|83p and ([5^ . This is because 
in our approach, the uncertainty is entirely due to sup- 
pressed higher order corrections. 

Finally, we note that our analysis does not account 
for the M^-dependence of the momentum-dependent, sub- 
leading contact interactions in the chiral EFT Hamilto- 
nian. Assuming validity of the naive dimensional analysis, 
such effects are beyond the accuracy of our present work. 
In addition, the strong correlations we observe for the M^- 
dependence of the various observables (as discussed in Sec- 
tion [7]) indicate that the relevant dynamics is largely gov- 
erned by the large S'-wave NN scattering lengths. Higher- 
order M^-dependent short-range terms are therefore ex- 
pected to play a lesser role. 



7 Correlations and the binding energy of the 
alpha particle 

Given the results in Section [SJ we are now in a position 
to draw conclusions concerning the individual energies E^ 
and the associated energy differences. The first interesting 
observation is that the energy differences AEy^, ^'^'^ 
e are, by themselves, extremely sensitive to changes in 
as could be expected. Such a conclusion follows from 
the unnaturally large coefficients in Eqs. (|55 )) - (j65[) . Notice 
that the fact that AE^^ , AE^^ and e are much smaller than 
the individual i?j does not, by itself, imply a strong fine- 
tuning. For example, the sensitivity of the Hoyle state 
excitation energy AE^ in Eq. (|73p is of a natural size, 
in spite of Z\i?™P — 7.65 MeV being almost an order of 
magnitude smaller than |£'i2^| = 92.16 MeV. 

While we find that AEy^ , AE^^ and e are by themselves 
extremely sensitive to variations in iV/^, we also observe 
the approximate relations 



h 



1.25, 



1.05, 



(84) 



which are satisfied for the central values of the individual 
terms in Eqs. at the level of a few percent. This 




Fig. 5. Sensitivities of £g, E12, -B12 and AE^ to changes in 
M^, as a function of K'^^ under independent variation of 
and over the range { — 1 ... 1}. The bands correspond to E12, 
E12, Eg and AE,. in clockwise order. 



suggests that AEf^, AE^^ and e cannot be independently 
varied (or fine-tuned) by changing the singlet and triplet 
NN scattering lengths (or, equivalently, by changing the 
strength of the short-range NN force in the ^5*0 and '^^i 
channels). Moreover, it is apparent from Eqs. (|55)) -([5 ^ 
that also the energies E^ of the individual states are as 
well strongly correlated in a similar manner. 

In order to quantify and illustrate the observed corre- 
lations, we show in Fig.[5]the changes in the sensitivities of 
iJg, Ei2, E12 and AE^ as a function of K"^^, when A^ and 

A^ are independently varied over a large range. The cor- 
relations for AEf^, AEi^ and e are shown in a similar way 
in Fig. [51 Within the statistical accuracy of our AFQMC 
results, we may conclude that the scenario of independent 
variations of the energy levels pertinent to the triple-alpha 
process under changes in the fundamental parameters is 
strongly disfavored. Given the prominent role of the '^He 
binding energy in the correlations, the observed behav- 
ior is strongly suggestive of the a-cluster structure of the 
^Be, ^^C and Hoyle states. Such correlations related to the 
production of carbon have indeed been speculated upon 
earlier ^MM- 



8 Reaction rate of the triple-alpha process 

We now turn our attention to the reaction rate of the 
triple- alpha process as given by Eq. and determine the 
range of variations in and a^,^ compatible with the for- 
mation of significant amounts of carbon and oxygen in our 
Universe, and thus with the existence of carbon-oxygen 
based life. We recall that the stellar modeling calcula- 
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Fig. 6. Sensitivities of AEf^, AE,j and e to changes in M^, as 
a function of K'^ under independent variation of ^4^ and At 
over the range { — 1 ... 1}. The bands correspond to AEf,, e and 
AEf^ in clockwise order. 



imply that a change in of as little as ~ 0.15% would 
suffice to render carbon-oxygen based life unlikely to exist. 
Stated differently, the "survivability band" corresponding 
to \dmg/mg\ < 0.15% would cover the whole of Fig. [T] It 
should be noted that in such a generic scenario, one can 
approximate 



de 



5M, 



1.5 



(87) 



MS 



which implies that the binding energy of ^He should be 
fine-tuned under variation of to its empirical value 
at the level of ~ 0.25% in order to fulfill the condition 
(5(e) I < 100 keV. Nevertheless, there clearly also exists a 
special value for the ratio of to Af. , given by 



AjAt ^ -1.5, 



(88) 



for which the dependence of ^i^/j, ^-E^ and e on be- 
comes vanishingly small (compared to the statistical un- 
certainties of the AFQMC calculation) , such that the fac- 
tor 

0.572(19) A^ + 0.933(15) Af - 0.064(6) < 1, (89) 



tions of Refs. [T51IM] suggest that sufficient abundances 
of both carbon and oxygen can be maintained within an 
envelope of ±100 keV around the empirical value of e = 
379.47(18) keV. For small variations |(5q;(,,-„/q;(,,-„| ^ 1 and 
Idm^/m^l <^ 1, the resulting change in e can be expressed 
as 



6is) 



de 



de 



MS 



MS 



de 



(85) 



where we recall that K'j^j 



4Q4+0 009 



, Thus, the 

condition \5{e)\ < 100 keV together with Eq. ([701) leads 
to the predicted tolerance |(SQ;(,„-,/a(,„J ~ 2.5% of carbon- 
oxygen based life to shifts in a^^. This result is compatible 
with the ~ 4% bound reported in Ref. [12] • For shifts in 



m^, we find 



0.572(19) ^ 
< 0.15%, 



0.933(15)^4-0.064(6) 



5m„ 



(86) 



using Eq. ([62]) , where we have neglected the relatively in- 
significant errors introduced by and The resulting 
constraints on the values of A^ and A^ compatible with 
the condition \S{e)\ < 100 keV are visualized in Fig. [7| 
The various shaded bands in Fig. [7| cover the values of A^ 
and Af consistent with carbon-oxygen based life, when 
is varied by 0.5%, 1% and 5%. 

In the most generic scenario, assuming that both of 
the dimensionless quantities A^ and A^ are ^ 0{1), and 
therefore 0.572(19)^, + 0.933(15) - 0(1), our results 



in Eq. (|86)) . In this case, we would conclude that the re- 
action rate of the triple-alpha process were completely in- 
sensitive to shifts in m^. As the realistic scenario is likely 
to be found somewhere in between these extreme cases, 
it becomes important to consider the available constraints 
on A^ and A^ before final conclusions are drawn. 

The quark mass dependence of the S'-wave NN scat- 
tering lengths has been analyzed within the framework of 
chiral EFT by several groups. The problem common to 
these calculations is the lack of knowledge about the m^- 
dependence of the NN contact interactions. Estimating 
the size of the corresponding LECs by means of dimen- 
sional analysis typically leads to large uncertainties for 
chiral extrapolations of the scattering lengths. For exam- 
ple, the NLO calculation of Refs. [401141] resulted in the 
values 

X9 =5±5, A'«^ = 1.1 ±0.9, (90) 

for the relevant JiT-factors. These are consistent with the 
NLO analysis of Ref. [12] , which yielded 



i^^ =2.4±3.0, i^^ =3.0 ±3.5, 



(91) 



based on a perturbative treatment of the OPE (see Ref. [2] 
for a related study). More recently, attempts have been 
made to combine chiral EFT with lattice QCD calcula- 
tions. In particular, the NPLQCD collaboration has de- 
termined the regions for the 5-wave scattering lengths 
consistent with their lattice results, — (0.63 ± 0.50) fm 
and at = (0.63 ± 0.74) fnr, obtained for M„ = 353.7 ± 
2.1 MeV [13]. By using these lattice data in conjunction 
with the assumptions of perturbativeness of the OPE po- 
tential in the ^Si-^Di channel and validity of the chiral 
expansion for NN scattering for > 350 MeV, Refs. [IS] 
H5] obtained results for and which are consistent 
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with the ones quoted in Eq. (|90|) . and in sHght disagree- 
ment with those of Eq. ((9T|) . 

Very recently, the analysis of the m^-dependence of 
NN observables was extended to N^LO in chiral EFT fIE\ . 
To overcome the difficulties due to the poorly known to^- 
dependence of the short-range NN interactions, the au- 
thors of Ref. [55] exploited the fact that the LECs ac- 
companying the NN contact interactions are saturated by 
exchanges of heavy mesons [47| . By means of a unitarized 
version of ChPT in combination with lattice QCD results, 
which describes the m^-dependence of the meson reso- 
nances saturating these LECs, the m^-dependence of the 
NN observables was analyzed at N^LO without relying on 
the chiral expansion of the NN contact interactions. The 
most up-to-date values are then given by 

Kl = 2.3tH Kl = 0.32t°;}^, (92) 

and 



^* = ~VW^ " -o.i8^°o;}°o> (93) 

which are not only consistent with the earlier determi- 
nations, but also in reasonably good agreement with the 
(parameter-free) LO chiral EFT calculation of Ref. [35], 
which is based on a novel, cutoff-independent approach. 
Interestingly, direct application of the central values in 
Eq. gives A^/ ~ —1.6, which leads to a strong can- 
cellation of the dependence of Zi-E/j, ^-Ef, and e on 
and , and hence also to a mild dependence on . The 
range of values given in Eq. (j93|) suggests that all contri- 
butions to the iiT-factors (notably including ) arc of 
^ 0(1), with the exceptions of K"^^ , and K^. 



9 Summary and conclusions 

We may summarize our findings by a brief discussion of 
Fig.[71 where we have superimposed the result correspond- 
ing to Eq. (j93p with the "survivability bands" from our 
AFQMC results in Eq. (|86|) . Given the current theoretical 
uncertainty in A^ and A^., our results remain compatible 
with a vanishing de/dM^, in other words with a complete 
lack of fine-tuning. Interestingly, Fig. [7] also indicates that 
the triple-alpha process is unlikely to be fine-tuned to a 
higher degree than ~ 0.8% under variation of m^. The cen- 
tral values of A^ and A^ from Eq. suggest that varia- 
tions in the light quark masses of up to 2 — 3% are unlikely 
to be catastrophic to the formation of life-essential car- 
bon and oxygen. A similar calculation of the tolerance for 
shifts in the EM fine-structure constant suggests that 
carbon-oxygen based life can withstand shifts of ~ 2.5% 
in a^j^. Beyond such relatively small changes in the funda- 
mental parameters, the anthropic principle appears neces- 
sary to explain the observed abundances of ^^C and ^^O. 



1| 




A 

s 

Fig. 7. "Survivability bands" for carbon-oxygen based life 
from Eq. (|86|l . due to 0.5% (broad outer band), 1% (medium 
band) and 5% (narrow inner band) changes in in terms of 
the input parameters Ag and At. The most up-to-date N^LO 
analysis of A^ and Af corresponding to Eq. (|93p is given by the 
data point with horizontal and vertical error bars. 

We also note that the fine-tuning in the fundamental pa- 
rameters is much more severe than the one in the energy 
difference e. 

Our ab initio lattice calculations account for all sources 
of quark mass dependence (explicit as well as implicit) in 
the LO nuclear Hamiltonian in chiral EFT, although we 
have not performed a strict LO analysis of the triple-alpha 
reaction rate. We have considered the potential impact of 
neglected higher-order terms on our results, in particular 
that of the 3NF which starts contributing at N^LO in 
the chiral expansion, and found that our conclusions are 
likely to be robust against such effects. Therefore, the most 
immediately useful extension of our work would be the 
incorporation of a more precise determination of A^ and 
A^ from future lattice QCD studies. 

As a longer-term objective, we may envision the in- 
clusion of dynamical photons in our AFQMC framework. 
Such a coupling of lattice QED to lattice chiral EFT may 
provide a more fundamental understanding of the sensi- 
tivity of the triple-alpha process to shifts in a^,^. 
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